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Abstract 



In this paper dividing factors Gl and Gt are constructed for the longitu- 
dinal and transverse responses of the relativistic Fermi gas in such a way that 
the reduced responses so obtained scale. These factors parallel another divid- 
ing factor studied previously, Hl, that yields a (different) reduced response 
which fulfills the Coulomb sum rule. Gl, Gt and Hl are all found to be only 
very weakly model-dependent, thus providing essentially universal dividing 
factors. To explore the residual degree of dependence which remains, the scal- 
ing and sum rule properties of several specific models have been considered. 
It is seen that the relativistic Fermi gas (by construction) and also typical 
shell-model reduced responses successfully scale and satisfy the Coulomb sum 
rule, as do experimental results at medium to high momentum transfers. On 
the other hand, it is observed that the quantum hadrodynamic model does so 
only if interaction effects become weaker with increasing momentum transfer, 
as predicted in the most recent versions of that model. 
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I. INTRODUCTION: REDUCED RESPONSE FUNCTIONS 



Motivated by our recent work on the Coulomb sum rule (CSR), in this paper we 
return to the related problem of ^-scaling. In Refs we devised a factor Hl that, upon 
being divided into the longitudinal response function for quasielastic electron scattering, Rl, 
produced a so-called reduced response, 



with convenient energy-weighted moments. Our approach started from the relativistic Fermi 
gas (RFG) model in which the moments are well understood. Specifically, the zeroth energy- 
weighted moment is the familiar CSR and becomes unity, that is, by construction is saturated 
exactly in the non-Pauli-blocked region where the momentum transfer q > 2kp, with hp the 
Fermi momentum. While the derivation of Hp provided in our previous studies was initially 
undertaken within the context of the RFG, we also showed that this dividing factor is only 
very weakly model-dependent, i.e., is essentially universal. 

A first motivation in the present work is to re-cast the ideas involved in studying the 
behavior of the quasielastic response in the region where the energy transfer uj is lower than 
its value at the quasielastic peak, namely, in the so-called scaling region. Below we start 
by summarizing the basic ideas behind the concept of y-scaling in which one attempts 
to find some function of q and u, here denoted G, which, when divided into the inclusive 
electron scattering cross section, yields yet another reduced response with special properties. 
Namely, for an appropriately chosen scaling variable y (a well-defined function of q and u; see 
below), this reduced response is a function only of q and y, denoted F(q, y), and scales. The 
latter means that for sufficiently large momentum transfers the function becomes universal, 
namely a function only of the scaling variable y: 



The arguments for choosing the dividing function and scaling variable may be presented 
from various points of view, always with the goal of removing the single-nucleon content 
from the nuclear responses in as model-independent a manner as possible while still retaining 
essential relativistic effects whenever feasible. In Sec. IIA we review the usual approach |3j 
based on the Plane- Wave Impulse Approximation (PWIA). This yields the standard form of 
y and F(q, y). Following this review, in Sec. IIB we re-cast our previous treatments of scaling 
H in terms of a related dimensionless variable if) which arises naturally when studying the 
RFG model. In the present work we show how if) is directly connected to a dimensionful 
variable t/rfgj which for the RFG model is the analog of the usual y-variable. Thus, scaling 
behavior can be examined in terms of y, ?/rfg or if). Each approach introduces specific 
functions by which the inclusive cross section or the individual longitudinal and transverse 
response functions, Rl and Rt, are to be divided to yield reduced responses F, Fl and Ft, 
respectively. In analogy with the previous treatment of the CSR, these dividing functions 
are denoted G, Gl and Gt, respectively; that is, for example 



r L = R L /H L , 



(1.1) 



F(q,y) 



F(y) = F(oo,y). 



(1.2) 




(1.3) 
(1.4) 
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In presenting results in Sec. Ill we show that, as was found to be the case for the dividing 
factor Hl, the dividing factors obtained in the scaling analysis are also only very weakly 
model-dependent. Indeed, in recent work medium- and high-energy data have been tested 
with both the usual y-scaling approach || and also the RFG-motivated approach || and 
seen to scale successfully (see also [7J for an analysis of all existing data). Additionally, it 
now appears that the experimental CSR is reasonably well saturated at high momentum 
transfers M. 

The second motivation for the present work is to put several models to the test, exam- 
ining both their CSR behavior and their ability to scale. Clearly the RFG model has both 
properties, by construction. We also examine two other models of the quasielastic response 
(see Sec. IIC), the hybrid model (HM) introduced in Ref. [|l] and the quantum hadrodynamic 
model (QHD) discussed in Refs. PJlQfl, to test how well they satisfy the CSR and are able 
to yield scaling behavior at high q. Each contains specific types of interaction effects that 
go beyond the strict RFG model and thus some violation of scaling and inability to yield 
a CSR of unity are to be expected. Our goal is to quantify the size of such effects and to 
explore whether or not the (successful) experimental behaviors can be used to constrain the 
models (see Sec. III). 

Finally, in Sec. IV we return to summarize our findings from this study. 



II. FORMALISM 

A. The reduced response function F and scaling 

We begin the discussion of the formalism with the most familiar chain of logic in which 
scaling is motivated within the context of the PWIA for (e,e'N) reactions HQ. Since in 
this approach one uses the integrals over (e, e'p) and (e, e'n) cross sections to approximate 
the inclusive (e, e') cross section, one may begin by performing the average of the electron- 
nucleon cross section over the azimuthal angle of the ejected nucleon: 

a eN (q,uj;p,£) = — f d(p N a eN (q,u;p,£;(p N ), (2.1) 

where p = \p\ is the missing momentum and £ the missing energy up to an offset of the 
(constant) separation energy E$ = + M < ^_ 1 — M\. Here is the nucleon mass, M\ 
is the target mass and M^_ 1 is the mass of the daughter nucleon when in its ground state 
(see Ref. [|IJ for the notation adopted in this work). Traditionally one uses some off-shell 
prescription for the electron- nucleon cross section, e.g., the "cc" prescriptions of De Forest 

n . 

Next one defines the scaling variable. The exact kinematics (i.e., no PWIA modeling is 
involved — only the imposition of energy-momentum conservation on the (e, e'N) process 
is required) permit one to say that the smallest value of the missing momentum attained in 
the so-called y-scaling region, the low-u; side of the quasielastic peak, occurs at p m in = —y, 
where || 
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W 2 +(m a _ 1 )'' 



m 



N 



(2.2) 



with 



W = \f(M A + u) -q 2 . (2.3) 



The energy transfer is, of course, then given as a function of q and y. In particular, it must 
lie in the range [l[ uo t < uo < q, where 



u t = E s + VWa-i + ™n) 2 + q 2 ~ (M°_ x + m N ) (2.4) 
is the threshold energy. The scaling variable vanishes when 



to = lu = E s + \J m 2 N + q 2 - m N , (2.5) 

which is roughly the position of the quasielastic peak, and hence the scaling region is char- 
acterized by having y negative. It should be stressed that no approximation is involved in 
using (q,y) as the inclusive scattering variables rather than (q, u), as they are related in a 
well-defined way through Eqs. ( [2.2| , [2.3| ). In discussions of quasielastic scattering the former 
usually proves to be more convenient, since the peak is found near y = 0, providing the best 
reference point for energy transfer at constant momentum transfer. 

The next step usually involves making two approximations. In the first, one assumes that 
protons and neutrons are distributed the same way in the nucleus, presumably a reasonable 
approximation for the N = Z nuclei considered in this work. Consequently one is assuming 
at this point that the spectral function of the nucleus does not carry any isospin index and 
that all such dependence can be incorporated in the overall eN cross section, 

E 

a eN (q,u;p,£) = — {Za ep (q,u;p, £) + Na en (q,u;p, £)} , (2.6) 

where the kinematic factor E^/q with E^ = ((<? + p) 2 + m 2 N ) 1 / 2 has been included, as in 
Ref. [1. 

Following this one attempts to remove this effective eN cross section from under the 
integrals over missing energy and missing momentum involved in going from coincidence to 
inclusive scattering. The usual argument is to assume that the most important contributions 
to the nuclear spectral function arise for the lowest values of (p, £) that can be reached for 
given values of q and y; in the scaling region these are £ = and p = —y. Thus, in this 
approximation, the function one hopes will scale as a function of y when q —>■ oo is given by 

°eN{q,y;p= -y,£ = 0) 

In the PWIA where the coincidence cross section factorizes into the electron-nucleon cross 
section and the nuclear spectral function S this quantity becomes 

Y 

F(q,y) = 2?r J p dph(q,y;p), (2.8) 
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where 

n(q,V',p)= I dSS(p,S). (2.9) 



o 



Here the conventional definition of the spectral function's normalization is taken to be 

roc _ 

dp d£S(p,£) = l. (2.10) 



o 



For the definition of the upper limits of integration in Eqs. Q2.8Q and ( |2.9| ) we refer the reader 
to Refs. [U and 0. Here we quote only the asymptotic results: 



\im Y(q,to) = oo (2.11) 

q^oo 



and 



q— >oo 



Jim (<?, y ; p) = y + p - xM-i + ~ V m a-i + 2/ 2 • ( 2 -l2) 



One sees that in the limit as q — > oo the result in Eq. (|2.12|) becomes independent of q and 
hence that the expression in Eq. (|2.8|) becomes a function only of y, namely, it scales. 



B. The reduced responses F^ FG and r^ FG 

Turning now to the RFG model, we have for its spectral function 

3A 



8nkl 



9{k F -p)8 S(p)-S iibG (p) 



where 



£ RFG (p) 



h\ -\- tn\; 



p 2 + m 2 N 



(2.13) 



(2.14) 



and where the usual convention is followed here of normalizing this spectral function to 
Z = N = A/2: 



dpi d£S RFG {p,S) = A/2 



(2.15) 



(note the different normalization from the one assumed above). For the physical responses 
discussed below, at the end of any calculation we, of course, take A — * oo in the strict 
RFG model. Defining the RFG scaling variable through the intercept of the support of the 
RFG spectral function given in Eq. (|2.13|) and the kinematical boundaries in the missing 
energy-missing momentum plane (see Ref. |T]|) we obtain 



2/rfg = m N ( = m N I AW 1 H k | , 



(2-16) 
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where, as in our past work, the dimensionless variables A = o;/2mjv, k = qjlm^ and 
t = k 2 — A 2 have been introduced. The RFG scaling variable in Eq. (|2.16| ) vanishes for 



Vl + 4k 2 - 1 , (2.17) 



namely at (uj — Es)/2m^, the value previously obtained shifted by the separation energy 
(usually a small effect, and one that can be incorporated as in our past work The 
different choice of scaling variable made here is motivated through the observation that, as 
we have chosen to describe heavy nuclei rather than few-body systems, we expect the main 
strength in the nuclear spectral function to be found not at the lowest values of p and E 
that are accessible, as is indeed the case for few-body nuclei, but rather for values of (p, 8) 
determined (largely) by the dynamics of nucleons in other than the ls-shell. 

Let us next close a logical loop and connect the RFG scaling variable defined above 
to the one originally proposed for the RFG || and referred to as ip. The latter, in the 
non-Pauli-blocked region, has the following equivalent forms: 



$ = -±= [20(A - A ) - 1] J(l + A) 2 + -(r - A) 2 - (1 + A) 



' [20(A-Ao)-1] Wl + --(1 + A) 



{ V r 

1 A-r 



(2.18) 



where £p — ep — 1 = yl+rjp — 1 and t]f = kpjvn^ are the dimensionless Fermi kinetic 
energy and momentum, respectively. Both ip and ( vanish for A = Ao and, by construction, 
for the RFG ip is confined to lie in the range — 1 < ip < 1. It is straightforward to see that 

$ = C ^ = (2-19) 

y/tF(l + VTT?) 

which implies that 

^ 2 = yfl + C 2 - 1. (2.20) 

The physical significance of ip is then immediately apparent: among the nucleons responding 
to an external probe one has the smallest kinetic energy and this is given by tp 2 (in units of 
the dimensionless Fermi kinetic energy £p). 

We may now proceed as in the previous subsection and divide by a factor that is pro- 
portional to the sum of the ep and en cross sections, weighted by Z and N, respectively. 
From our previous work on the RFG model || we know that the azimuthal-angle-averaged 
eiV cross section (compare Eq. (|2.1|)) is proportional to vl/l + vpfr, where vl and are 
the usual longitudinal and transverse Rosenbluth kinematical factors and, in terms of the 
electric Ge v ^ and magnetic Gjup.n Sach's form factors, we have 
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K 2 



T 



f L (K, A; X ) = - Z(G 2 Ep + W 2pX 2 ) + N(G 2 En + W 2nX 2 ) (2.21) 



/ t (k, A; X ) = Z(2rG 2 Mp + W 2pX 2 ) + N(2rG 2 Mn + W 2nX 2 ) (2.22) 



with 



W 2p , n (r) = — [G 2 Eptn (r) + rG^ n (r)j . (2.23) 

Here x = ??sin#, where 77 = p/m^ and # is the angle between p and q. Note that the 
5-function in Eq. ( |2.13|) requires that S = S RFG as defined in Eq. fl2.14|) . 



We could proceed to follow exactly the procedures outlined in the previous subsection, 
working from the unseparated inclusive cross section towards a reduced response which, if 
successful, would scale as q —>■ 00. Instead, in this work for the most part we work directly 
with the separated longitudinal and transverse responses, Rl and Rt, respectively, since 
(1) we are most interested in model-to-model comparisons and the same procedures may be 
followed in each case (i.e., focusing on L or T responses directly), (2) a few cases exist where 
L/T separations have been performed experimentally, and (3) we wish to draw comparisons 
with studies of the Coulomb sum rule where only the L response is relevant. In fact, in most 
of the discussions to follow we shall limit our attention to the longitudinal channel and only 
give a few results for the transverse case. 

We now seek reduced responses denoted Fl^t(k,iP) that scale. These are to be obtained 
from the inclusive response functions Rl,t{k,X) by dividing through by specific functions 
denoted Gl,t(k, A) (discussed below): 

F LjT (k, 4) = R LiT (k, X)/G L)T (k, A). (2.24) 

If the dividing functions are chosen appropriately, then as above the reduced responses 
defined in Eq. ( [2.24j ) will scale, namely, become functions only of a single scaling variable 
such as ip defined above when k — > 00, 

F l , t (k, ™ F LiT (iP) = F LiT (oo, V). (2.25) 

To obtain the dividing functions we proceed from the spectral function for the RFG, 
using the general PWIA expressions given above. Let us first focus on the longitudinal 
response. With Eq. ( |2.13|) we have 



3A 



/•■in 
d£S RFG (p, S) = ^0(k F - p)0(p + y RFG ) (2.26) 

F 

and therefore, from Eq. ( |2.8| ), together with a factor Z = N = A/2 to account for the 
different normalizations of S and S R , we have 

kp 

Fl FG (q,y RFG ) = Fl FG (y RFG ) = ^ J p dp = ^(k 2 F - y 2 RFG ) . (2.27) 

-wrfg 

By exploiting Eq. ( |2.2(J| ) it is then an easy matter to obtain 
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k 2 F - 2/rfg = mj^ F (l - ^ 2 )[2 + + ij 2 )] 



(2.28) 



and hence the scaling function of the RFG will read 

3£f 



jpRFG 



2m N rf F 



(i-W-f) 



(2.29) 



We now recall from previous work [||] that the RFG longitudinal response function is 
given by (see Fig. [I]) 



R 



RFG, 



K, A) 



36 



where the nucleonic terms Ui, p ,n are the following: 



{ZU Lp + NU Ln ](l-?p 2 )6(l-i; 2 ), 



(2.30) 



U. 



Lp,n 



K 
T 



G 2 Ep , n (r) + W 2p>n (r)A 



Furthermore A, which, in the non-Pauli-blocked domain, reads 



K 



-(e^ + e FV /l + C 2 + l + C 2 



+ A (ep + Jl + C 2 + A 



-1 + r) 



(2.31) 



(2.32) 



represents the transverse component (with respect to q) of the momentum of the struck 
nucleon [0. It follows that the required RFG dividing function is given by 



G l {k,X) 



zu, 



Lp 



NU. 



Ln 



2k[1+? f (1 + *' 2 )/2] 



2^ 



NU 



Ln) 



0(6 



(2.33) 
(2.34) 



telling us that the scaling function of the RFG arises from Rl not only by pulling out the non- 
scaling single-nucleon factor involving Ul p /k and Ui n / 'k, but also contains a small correction 
for the medium dependence in the problem (the factor in the denominator involving ip). 
Typically £p is very small, attaining values as large as 0.04 only in the heaviest of nuclei, and 
thus both the correction in the denominator containing £p and the medium-dependent effects 
in the numerator embodied in A provide only few percent corrections at high momentum 
transfers, as long as ip is not permitted to become large. Indeed, in the RFG its magnitude 
is bounded by unity and so this always obtains at high q; if the ideas here are carried over 
to other models where large excursions away from the quasielastic peak are permitted, then 
one should check the actual size of these medium-dependent corrections. 

Dividing the longitudinal response function by Gl then yields the reduced response F L 
which, at least for the RFG, scales with ip, ( or ?/rfg- Note that we can use Eq. (|2.20|) to 
write 



i + -Mi + ^ 2 ) 



2 {€F 



(2.35) 



giving us some insight into the physical significance of the denominator in Eq. ( |2.34j) : namely, 
it represents an average energy of the nucleons responding to the external probe. 
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Proceeding in exactly the same way for the transverse channel it is straightforward to 
see that the dividing factor required in discussions of scaling in that case is 

Gt(k, A) = ZUtp I NUt ^ , , (2.36) 

= ^-(ZU Tp + NU Tn ) + <D{£ F ), (2.37) 

where (compare Eq. (|2.22|) ) 

U Tp , n = 2rG 2 Mp:n (r) + W 2p , n (r)A. (2.38) 
The unseparated responses may then be analyzed using 

, s d 2 cr/dQ e du 

<tm[vlG l (k,\)+vtG t (k,\)Y 

where as usual one would have for the cross section 

d 2 a/dVL e duj = ct m [v l R l (k, A) + v t R t (k, A)] (2.40) 

with crjvf the Mott cross section and Vl 7 t the usual Rosenbluth kinematical factors. Naturally, 
in the RFG model we have F RFG = F^ FG = F^ , the result being given in Eq. (|2.29| ). 

The developments presented here for scaling have been motivated by our previous dis- 
cussions of the Coulomb sum rule and various energy-weighted moments of another reduced 
response denoted Tl(k, A) (]],§]. In those studies the longitudinal response Rl(k,X) was 
divided by a function H l (k, A) to yield 

t l (k,iI>) = R l {k,\)/H l (k,\), (2.41) 

where the ia th moment of the longitudinal response of the nucleus is given by 



d\X n r L (K,\). (2.42) 



In particular, the n = moment, 5^, is the Coulomb sum rule. In the case of the RFG our 
previous work showed that 

ZU Lp + NU Ln 

H l (k,X) = , (2.43) 

Jl 

where Ul PjTI are given above and where 

J, = f^ (2.44) 
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Thus, by dividing the charge response of the RFG by Eq. ( |2.43| ) (see Eq. ( |2 . 30|) ) one obtains 
the following reduced longitudinal response 



„RFG , 



(2.46) 



which, by construction, fulfills the CSR in the non-Pauli-blocked domain, as can easily be 
verified. Observe also that while Fj^ FG attains its maximum for ip = 0, i. e. for A = Aq, this 
is nearly, but not exactly, true for r^ FG only for small values of r] F . Indeed using the simple 
expression 



if, 



1 

Vf 



A(A + 1) 



K 



(2.47) 



which approximates Eq. (|2.18|) quite faithfully for very large values of k except on the borders 
of the response region, it turns out that the maximum of r^ FG occurs at 



A 



\ 



1 + An 2 1 



1 Vf 
V5 



A + 0{rf F ) 



(2.48) 



Before leaving this discussion of the RFG model reduced responses, for completeness 
let us mention another version of the reduced longitudinal response that can be obtained 
in terms of the variable ( rather than ip as in Q2.41| ). For this purpose one has to divide 
Eq. ( CT) by Eq. (gg3| ) with J L replaced by 



_ 2nrf F d( 



3 d\ sinh 1 rj F 
e F - 



(2.49) 



Vf 



The result is 



„RFG' 



K,X) 



O r T^l 9 WF - C ) oT", 

2m N r] F smh r] F OA 



(2.50) 



e F 



Vf 



which again fulfills the CSR. Indeed 

1 



^(0),RFG 



dXrf FG 'i 



k, A) 



Vf ( sinh r] F s 
e F 



\-i v r ^ ( €F ~ V / l+^) = L ( 2 - 51 ) 
mn rip \ J-rsp \ J 



Vf 



C. Scaling and sum rules in the HM and QHD models 

An extension of the RFG to account for the binding of the nucleons inside the nucleus, 
thus curing a flaw of the RFG related to its negative separation energy, is represented by the 
hybrid model studied previously in Ref . . The HM has continuum states which are simply 
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plane waves, as in the RFG model, but has bound states described by shell-model wave 
functions obtained by solving the Schrodinger equation with some choice of potential well. 
For simplicity in studying the A — > oo limit, in M we chose to consider harmonic oscillator 
bound-state wave functions. The longitudinal response R^ M is also shown in Fig. p]. 

In the HM scaling variables which incorporate the shift result from using Eq. ( |2.16| ) 
for (and a corresponding dimensionful variable y') or Eq. ( |2.18|) for ip' by making the 
replacements A — > A' and r — > r' = k 2 — A /2 , where 

A' = A - A shift , (2.52) 

being 

\bi&=-^—(T F + E s ) (2.53) 



2m 



N 



and Tp = y k F + m N — rriN = win^f the Fermi kinetic energy. The HM turns out to have the 
width of its reduced response (its variance, see Ref. |3J) identical to a RFG computed with 
a Fermi momentum somewhat larger than the usual one (237 MeV/c for the HM, versus the 
230 MeV/c value for kp used here for the RFG to correspond to nuclei near 40 Ca or 56 Fe). 
In conformity, ip' should be computed with kp=237 MeV/c as well. As seen in Fig. p], the 
shift arising from replacing A with Ashift according to Eq. ( j2.53| ) is apparent. 



Scaling may then be examined for the HM by computing 

Fl M (n, 4) ee R™(k, X)/G l (k, A) (2.54) 

and the various energy-weighted moments of the longitudinal response, including the zeroth 
moment or CSR, may be computed using 

r™M) = RT(KA)/H L (K,\). (2.55) 

Note that we have used the same dividing factors Gp and Hp that were developed from our 
discussions of the RFG in the previous subsection. In [JTJ] we showed that H L is universal 
in that the form found for the HM essentially coincides with the one obtained above. In 
the next section we present results where the corresponding universality of Gp is tested 
numerically. Naturally these explorations also involve displaying the reduced responses 
given in Eqs. ( |2.54|) and (|2.55| ) as functions of ip' as well as ijj. 

Finally we turn to an examination of the QHD model @. In this model protons and 
neutrons in the nucleus are described by Dirac spinors and move in strong Lorentz scalar 
and vector mean fields. These in turn arise self-consistently from the exchange of a and 
uj mesons between the same nucleons on which they act. The scalar field dresses the bare 
mass of the nucleon, considerably lowering its value; the vector field uniformly shifts the 
fermion spectrum. As a consequence the QHD charge response of nuclear matter in Hartree 
approximation is unaffected by the vector field, while it turns out to be quite sensitive 



to the effective mass m* N induced by the scalar field, as shown in Ref. fLOfl . This is, of 
course, true in the simple approximation of constant relativistic mean fields. An improved 
description allows for an energy dependence of the latter which helps to account for the 
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data of proton-nucleus elastic scattering. We shall however ignore these refinements in the 
present analysis. 

In Fig. in order to directly compare with the results of Ref. | ID| , we display the 



longitudinal response of nuclear matter, again using kp = 230 MeV/c, in the Hartree version 
of QHD for q = 0.55 and 1.14 GeV/c, respectively. In addition to the RFG and HM results 
discussed above, here we also show curves corresponding to m* N = 0.8 and m* N = 
0.68 rriN, namely, those displayed in Ref. []TU| for the case of constant relativistic mean fields. 
From the figure the sensitivity of the charge response to m* N is fully apparent. Clearly the 
more m* N deviates from the more the longitudinal QHD response is hardened (shifted 
to higher values of u). The most recent version of QHD suggests that as q increases the 
effective value of m* N should tend towards and thus that, for intermediate-energy studies 
(q ~ 500 MeV/c), m^/mjy may be somewhere between 0.68 and 0.8, while for high q, a 
value between 0.8 and 1.0 is likely to be preferred. 

As above for the HM model, scaling will be examined in the QHD model by computing 

F l qhd (k, if)) = RT B (k, X)/G l (k, A) (2.56) 

and likewise the various energy-weighted moments of the longitudinal response will be com- 
puted using 

r? HD (^) = itg HD (K, \)/H l (k, A), (2.57) 

where R^^^k, ip) is given below in Eq. (|2.63|) . Our contention is that the dividing factors are 
(at least to a very good level of approximation; see Sec. Ill) universal and accordingly here 
we have continued to use the same dividing factors Gl and Hl that were developed from our 
discussions of the RFG in the previous subsection. In discussing the QHD results in the next 
section we shall display the reduced responses in Eqs. ( |2.56| ) and (|2.57| ) as functions both of 
if) and also if)*, namely, the RFG scaling variable given in Eq. fl2.18|) with mjy replaced by 
m* N . We shall also briefly explore what happens when the mjy — > fn* N replacement is made 
in a particular way to obtain modified dividing factors, G* L and if£, defined below. 

Before proceeding to reduce the QHD response a delicate point arises in connection with 
the question: does the concept of effective mass also apply to the single-nucleon content 



of the longitudinal response? The answer is yes and no. In fact the approach of Ref. [IH 
requires the Dirac F\ and Pauli Fi/m^ form factors in the medium to be identical to those in 
free space (or, in other words, the nucleon current to be unaltered). Thus they are expressed 
in terms of the Sach's electric and magnetic form factors through the standard relations 



" lp,n 

and 



» = — |— [G Ep>n {r) + rG Mp>n {r)} (2.58a) 

1 + T 



F 2p , n (r) = — ^- [G Mp , n {r) - G Ep . n {r)} , (2.58b) 
1 + r 

with Ge p ,u and Gup,n fitting the results of electron-nucleon elastic scattering in free space. 
Then, to keep Fi PiU and F 2p , n /m N unaltered, the Sach's form factors in the medium should 
be given by 
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G* Epjn (q,u;) = F lp , n (r) - T *^-F 2p>n (r), (2.59a) 



G* Mp , n (q,uj) = F lp , n (r) + -^F 2p ,„(r), (2.59b) 

m N 

where r* = {m^ /m* N ) 2 T. Hence they will differ from the corresponding quantities in free 
space. Indeed it is the m* N appearing in Eq. ( [2.59b|) that causes a strong reduction of the 
proton's convection current and, as a consequence, leads a considerable quenching of the 
magnetic moment of odd-proton-nuclei. To restore the agreement between the theoretical 
predictions of QHD and experiment a substantial back-flow current contribution has then 
to be invoked [|T2| . 



For the quasielastic charge response we face a similar situation. In fact, as discussed 
above, we reduce the QHD longitudinal response by using as dividing factors Gl and Hl 
defined in Eqs. (|2.34j) and (|2.43|) , respectively, with the bare nucleon mass and, hence, with 
the free Ge p ,h and Gm p ,u- An alternative might be to use G* Epn and G* Mpn given above, 
together with 

W; p , n (q, w) = y^-t {Gt,n + r*Gtr p>n ) , (2.60) 

to define £/£ using an expression that is analogous to Eq. (|2.23|) . This would lead to dividing 
factors we denote G* L and and, by making the replacement — > m* N in Eq. ( |2.1£ ), to 



a corresponding scaling variable ip*. Accordingly we would have different reduced responses 
(again labeled with an asterisk): 

F L QHD >,^) = R2 nD (K,X)/Gl(K,X) (2.61) 

and 

r^M) = R^ D (k,X)/HI(k,X). (2.62) 
Clearly, since the QHD longitudinal response function is given by (compare with Eq. ( |2.30| )) 

i?? HD ( K , A) = -J®[ZU^ + NUU(1 - r 2 W(l ~ r 2 ) , (2.63) 

the reduced responses i^Q HD * anc i r Q HD * w j}} behave as functions of if)* exactly as i^- FG and 
r RFG bg^avg as functions of ip. In the following section we return briefly to examine the 
consequences of these alternative definitions of the dividing factors. 

III. RESULTS 

We have thus reached the point at which the longitudinal reduced response functions Fl 
(for scaling discussions) and (for sum-rule discussions) are to be obtained by dividing Rl 
by two different factors, Gl and Hl, respectively. The first is given by Eq. Q2.34Q , while the 
second is given in Eq. Q2.43Q and they are related according to 
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G ^- A) ^(i|)lx 1+ ^(l 1 + ^)/2 g ^ A> (31) 



1 /k\ ( 1 + 2A 



2 \tJ \ 1 + A 



A) +0(60- (3.2) 



Naturally, the analogous dividing function Gt in Eq. ( [2.371) should be used when treat- 
ing the separated transverse response or the combination (Tm[vlGl(k, A) + VtGt(k,X)) in 
Eq. ( |2.39| ) may be used when dealing with an unseparated cross section. To assess the 
universality of using the various approaches discussed in the preceeding section, we show 
in Figs. |2| and ^ two mixed representations of the usual y-scaling results and the RFG- 
motivated results. In Fig. || the total RFG cross section is displayed divided by the De 



Forest ccl off-shell result [T1J - this is what is usually set as the dividing factor in pre- 
senting unseparated data in the form F(q,y) versus y, as in Ref. ||. Here we have taken 
Eg = 8 MeV as being typical of nuclei across the periodic table. Clearly once the momentum 
transfer is large enough (say above about 1 GeV/c) it is irrelevant whether the ccl form is 
used or the RFG form is employed: the dividing factor is effectively universal. In Fig. |3| 
the usual reduced response for the RFG (i.e., with the dividing factor given in Eq. ( |2.39| )) 
is displayed, however as a function of y rather than if), in which it of course scales perfectly. 
Again, at all but the lowest momentum transfer the scaling is excellent in y as well as if). 
These and other more extensive studies show that at high momentum transfers the choice 
of scaling variable is largely irrelevant, as long as one of the relativistic forms discussed in 
this work is employed, and that the on- and off-shell prescriptions for the dividing factors 
yield essentially model-independent results for scaling (as they did for the CSR). 

By exploiting the essentially universal dividing factors Hl and Gl, we are thus in a 
position to test whether or not a given model fulfills the CSR and scales. Here we examine 
the longitudinal response functions of the HM |]J and QHD Pl,|l0l| models in concert with 



the RFG results obtained previously. 

The response function ir^ M of the hybrid model has been computed in Ref. and, when 
divided by H L , shown to satisfy the CSR. There, as already mentioned, it was seen that 
in setting up if)' the effective Fermi momentum 237 MeV/c should be taken: the associated 
shift in A (see Eq. (|2.53| )) turns out then to be A s hift = 0.02. Fixing Gl similarly, in Fig. |] 



we display F]^ M for 40 Ca versus if) and if)' for four different values of q (the RFG result is 
also shown for reference) and then display it in Fig. |5] versus q for if) = —0.1 and —0.5, as 
well as for if)' = —0.1 and —0.5. As seen in Fig. [|, the HM scales either with if) or with if)' 
as q becomes large; indeed, only the q = 500 MeV/c plot versus if) shows any appreciable 
violation of scaling. Clearly (by construction) the scaling versus if)' is excellent. It should 
also be noted that the final scaling result in the upper panel in Fig. |4] lies a little to the 
right of if) — 0, while that in the lower panel lies a little to the left of if)' = 0, suggesting 
that the value for A s hift in Eq. ( |2.53| ) used here is somewhat too large and has led to an 
"overshooting" of the shift. Figure |5| shows the evolution of scaling with q and confirms 
that the asymptotic behavior has essentially been reached by about 1 GeV/c. We also see 
that the scaling sets in sooner when \if)'\ is smaller, i.e., when one is closer to the quasielastic 
peak position. 

In Figs. |] and |7| we display i^ HD for 40 Ca versus if> and if)* for four different values of 
q (the RFG result is also shown for reference) and then display it in Fig. 18] versus q for 
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ip = —0.1 and —0.5, as well as for ip* = —0.1 and —0.5. Here results are shown both for 
m* N = 0.68 m N and 0.8 m N . The results in Figs. || and [7| show that i^ HD does not scale 
versus ip when the effective mass is constant and differs from m^. As q continues to grow 
beyond the range of values shown in the figures, the results continue to shift to higher uj and 
never coalesce into a universal curve (see also Fig. |8]). When plotted versus ip* the behavior, 
while better, still does not scale. This is contrast with the RFG and HM results displayed 
above and, importantly, is not what is seen experimentally where the world data do appear 
to scale in ip |7[ (see also ||). Comparing Figs. || and [7], it is clear that scaling is better 
when m* N is closer to m^, indeed becoming perfect when m* N — > mjv, as expected, since then 
pQ d _^ pRFG rj^g £ ac j. ^at experimentally the scaling is observed to occur successfully 
for q greater than about 1 GeV/c suggests that m* N /m N should not deviate appreciably 
from unity for such kinematics. 

In Fig. H the CSR associated with the QHD model is displayed for 40 Ca [kp = 230 
MeV/c) and for a set of values of m* N ranging from 0.5mjv to the bare nucleon mass. Here 
we observe an effective mass dependence mostly due to the factor m* N in front of the CSR 
integral 



•(°)>Q HD = / dujr^{q,u) 
o 



4 \m N ) ( F J V Ul ' 



Am N J U L dip 



but also, to a less extent, to the Jacobian dip j 'dip* . In fact, using Eq. ( |2.47| ), the latter can 
be cast in the form 

^ i _ l-m* N /m N ^ ^ 

9^* Jl + AK(K + T] F lP) ' 



which is less than unity for finite k. These two factors, which are both clearly related 
to many-body aspects of the nuclear response, are however partially counteracted by the 
"single- nucleon term" U^/Ul which is also m* N dependent and greater than unity. For large 
values of k, Eq. (|3.4j) goes to one, whereas UI/Ul continues to grow rapidly, increasingly so 
for smaller values of m* N . This then becomes the dominant factor in explaining the behavior 
of the curves in Fig. |9] for q, say, larger than 1 GeV/c. 

As discussed above, the expectation is that m* N must evolve with q towards m^. Thus, 
while at q — 500 MeV/c a value of m* N /rriN — 0.7 may be acceptable (implying that the 
resulting CSR will be about 5-10% below unity), as q reaches the higher values between 1 
and 2 GeV/c, values of m* N /m^ nearer unity will be preferable from the point of view of 
scaling (and also suggested by the full version of QHD, see Ref. [0) and will therefore yield 
a CSR that still remains close to unity. 

For completeness, using the nomenclature introduced in Ref. |] , we display in Fig. [LCj the 

CSR (H^ )), the energy-weighted sum rule (SW/S^) and the variance a = a/H( 2 ) — (HW) 2 
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of 40 Ca according to the HM, QHD (with m* N = 0.68 m N and 0.8 m N ) and RFG models. As 
expected the responses of the HM and QHD models are considerably hardened with respect 
to that of the RFG; in addition the width of the QHD response is substantially wider than 
for the other two models. 

In concluding this section we shortly discuss the scaling variables. In this paper four 
scaling variables have been introduced: the canonical y, the ip of the RFG (together with 
its related variations, ( and ?/rfg), the ip' of the HM and the ip* of the QHD. How are they 
interrelated? To answer this question in Fig. |TT| we display their behavior as functions of uj 
for q = 1 GeV/c. The following features emerge from the figure: 

a) The behavior of the canonical y variable versus the mass number A is very rapid at 
small A and appears by A = 20 to have almost reached its asymptotic A = oo value. 

b) In accord with the findings of for large A, y is roughly given by Eq. ( |2.47| ) with 
A shifted downward by an amount grosso modo corresponding to Eg- Hence in the A = oo 
limit one does not recover the RFG y- variable from y, but a slightly shifted result. 

c) The variable ip' is also shifted downward with respect to ip, by an amount set by A s hift 
in Eq. ( |2.53| ), although we conjecture (see above) that this shift is somewhat too large. 

d) Finally, note the behavior of ip* with uj: it turns out to be shifted towards higher 
energies with respect to the other scaling variables. This goes in parallel with the energy- 
weighted sum rule of the QHD model whose reduced response is substantially hardened in 
comparison with other models here considered for q >500 MeV/c. 



IV. FINAL COMMENTS 

Our motivations in this work have been to study the Coulomb sum rule and scaling in 
a unified way, and to explore various models for the inclusive electromagnetic response of 
nuclei in the quasielastic region to see whether or not they simultaneously satisfy the CSR 
and scale. We have demonstrated how our previous work on the CSR may be generalized 
to include treatments of scaling: in both cases the response functions or the cross section 
are divided by specific functions, H L for the CSR and Gl, Gt or <Jm{vlGl + vtGt), as is 
appropriate, for scaling, to yield reduced responses with the desired properties, namely, a 
CSR or scaling behavior. We have studied the model-dependence of these dividing functions 
within the context of two specific models for the single-nucleon content in the problem and 
found it to be small, as long as we restrict our attention to high momentum transfers and to 
the vicinity of the quasielastic peak. Essentially, for such conditions, we have obtained nearly 
universal dividing functions with which any model can be tested and with which experimental 
data can be reduced. Typically we find corrections at high q which are characterized by 
If = (^f/^at) 2 , namely, reaching only ~8% for the heaviest nuclei. 

In the course of this study we have introduced and inter-related several types of scaling 
variables (y, if), (, ?/rfg> V*0- These variables are all closely related: they differ at most 
by small shifts introduced by the various models to account (at some level) for interaction 
effects. On the one hand, when experimental data are reduced using our dividing functions, 
they are seen to yield a CSR at high-g and to scale when plotted versus any of the scaling 
variables listed above. On the other hand, the models considered in this work yield various 
results. The relativistic Fermi gas (RFG) model and the hybrid model (HM) both saturate 
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the CSR at high-g. In contrast, the quantum hadrodynamical (QHD) model does so only if 
the effective value of m* N /m^ evolves with increasing q towards unity, as suggested by the 
latest version of the model. With regard to scaling, the RFG model (by construction) scales 
with ip (and of course with ( or 2/rfgj since they are intrinsically related). It also scales 
quite well versus y. The HM scales very well with ip' (by construction), but also quite well 
with ip. In contrast, the QHD model does not scale with any of the above variables if the 
ratio m* N /mN is constant and differs significantly from unity. One must conclude that the 
successful scaling behavior seen experimentally indicates that m* N must approach as g 
increases beyond about 1 GeV/c. 

Finally, we note that an exact CSR and exact scaling behavior should not be expected 
to occur. We certainly believe that interaction effects beyond those of the mean field (whose 
effects are incorporated at least to some extent in the models studied here) can play a role 
even at momentum transfers as high as 1 GeV/c. Furthermore, we expect that two-body 
meson exchange current effects can also play a role, for instance, making the longitudinal and 
transverse reduced responses scale to different functions. Both of these classes of corrections 
are currently being re-investigated. 

ACKNOWLEDGMENTS 

This work was supported in part by funds provided by the U.S. Department of Energy 
(D.O.E.) under cooperative research agreement #DF-FC02-94ER40818 and by the INFN- 
MIT "Bruno Rossi" Exchange Program. 



17 



REFERENCES 



[1] R. Cenni, T.W. Donnelly and A. Molinari, Phys. Rev. C56, (1997) 276. 
[2] P. Amore, R. Cenni, T.W. Donnelly and A. Molinari, Nucl. Phys. A615, (1997) 353. 
[3] D.B. Day, J.S. McCarthy, T.W.Donnelly and I. Sick, Annu. Rev. Nucl. Part. Sci. 40, 
(1990) 357. 

[4] E. Pace and G. Salme, Phys. Lett. B110, (1982) 411; C. Ciofi degli Atti, E. Pace and 

G. Salme, Phys. Lett. B127, (1983) 303. 
[5] W.M. Alberico, A Molinari, T.W. Donnelly, E.L. Kronenberg and JW. Van Orden, 

Phys. Rev. C38, (1988) 1801. 
[6] C.F. Williamson et al, Phys. Rev. C56, (1997) 3152. 
[7] T.W. Donnelly and I. Sick, (to be published). 

[8] J. Jourdan, Phys. Lett. B353, (1995) 189; Nucl. Phys. A603, (1996) 117. 

[9] B.D. Serot and J.D. Walecka, Advances in Nuclear Physics 16, (1986) 1. 
[10] H. Kim, C.J. Horowitz and M.R. Frank, Phys. Rev. C51, (1995) 792. 
[11] T. De Forest, Jr., Nucl. Phys. A392, (1983) 232. 
[12] J. A. McNeil et al, Phys. Rev. C34, (1986) 746. 

[13] M.B. Barbara, A. De Pace, T.W. Donnelly and A. Molinari, Nucl. Phys. A569, (1994) 
701. 



18 




FIG. 1. The longitudinal response function Rl is displayed versus to for q = 0.55 GeV/c in 
panel (a) and q = 1.14 GeV/c in panel (b). The four curves shown are the following: dashed - 
RFG, solid — HM, dot-dashed — QHD model (m^=0.8 m N ) and dotted — (m^=0.68 m N ). Here 
and in all of the following figures the nucleus chosen is 40 Ca. For the RFG and QHD models we 
always assume /ci?=230 MeV/c, while the HM parameters are discussed in the text. 
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FIG. 2. Scaling is shown as a function of ip using the RFG cross section at 10° scattering 
angle divided by the off-shell ccl form with separation energy E$ = 8 MeV for a wide range of 
momentum transfers. 
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FIG. 3. The scaling function F RFG is shown for the same conditions employed in Fig. [2] as a 
function of the y-scaling variable defined in Eq. (|2.2| ), rather than as a function of ip (see Eq. (|2.18[) ) 
in which it scales exactly, by construction. The separation energy occurring in the definition of y 
has been set to 8 MeV for these results. 
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FIG. 4. The reduced response F^ M (the hybrid model scaling function) is shown as a function 
of ij) in panel (a) and tp' in panel (b) for 4 values of q (fine dotted — 0.5, dot-dashed — 1.0, heavy 
solid — 1.5 and solid — 2.0 GeV/c). The RFG result, which scales exactly as a function of ip, is 
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FIG. 5. The reduced response is shown as a function of g at tp = —0.1 and —0.5 in panel 
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and solid — 2.0 GeV/c). The RFG result, which scales exactly as a function of ip, is also shown 
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FIG. 7. As in Fig. |(| except now with m* N = 0.68 mjv 
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FIG. 9. Coulomb sum rule h(°)>Q hd of the QHD model is shown as a function of q for different 
values of the effective mass m* N , ranging from 1 to 0.5 tun in steps of 0.1 tun- 
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FIG. 10. The Coulomb sum rule (a), energy- weighted sum rule (b) and variance (c) are shown 
as functions of q for the three models employed in this work: dashed — the RFG, solid — the HM 
and the QHD model with m* N =0.8 mjy (dot-dashed) and 0.68 mjy (dotted). 



28 




200 300 400 500 600 700 



cj [MeV] 



FIG. 11. The various scaling variables employed in this work are shown as functions of uo at 
q = 1 GeV/c: k F ip (dashed), k F (HM)ip' with k F (HM) = 237 MeV/c (dot-dashed), k F tjj* (dotted) 
and y (solid). The three solid curves correspond to A = 2, 20 and 200, respectively, starting from 
below. 
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